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We show by general arguments that networks whose density of states contains few highly degenerate eigenval- 
ues result in inefficient performances of continuous-time quantum walks (CTQW) over these networks, while 
systems whose eigenvalues all have the same degeneracy lead to very efficient transport. We exemplify our 
results by considering CTQW and, for comparison, its classical counterpart, continuous-time random walks, 
over simple structures, whose eigenvalues and eigenstates can be calculated analytically. Extensions to more 
complicated, hyper-branched networks are discussed. 
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I. INTRODUCTION 

The transfer of information in quantum systems has at- 
tracted a lot of attention in recent years, especially in the con- 
text of quantum computing 1 1]. In analogy to classical random 
walks, which are used as algorithmic tool in "classical" com- 
puting, two version of quantum walks have been introduced: 
discrete-time quantum (random) walks, with an additional 
internal "coin" degree of freedom [2], and continuous-time 
quantum walks (CTQW), where the analogy to continuous- 
time random walks (CTRW) lies in identifying the classical 
transfer matrix with the quantum mechanical Hamiltonian [3]. 
Recently, it has been shown how these two version are related 
fl. 

CTQW are formally equivalent to the tight-binding model 
in solid-state physics pj or the Hiickel/LCMO model in phys- 
ical chemistry [6] and therefore can be applied to study trans- 
port processes in various types of different systems, like spin 
chains ff] or ultra-cold Rydberg gases |8]. What matters is 
that the constituting elements (spins, atoms, molecules, etc.) 
are of the same type, in the simplest cases they resemble two- 
level systems. 

In general, the quantum dynamics is very different from the 
coiTesponding classical dynamics. The unitary time evolution 
of a quantum system leads to characteristic dynamical phe- 
nomena, e.g., quantum revivals and quantum carpets for the 
particle in a box |9] or for CTQW on a ring 110(1 . or, for disor- 
dered systems, to localization [ 1 1.1 . 

For linear (ordered) systems, the quantum transport effi- 
ciency of CTQW has been proven to overcome the classical 
efficiency, a result which translates also to more complicated 
systems like decision or Cayley trees ifl^l . However, the effi- 
ciencyof the transport strongly depends on the initial condi- 
tion 1 13|], i.e., different initial conditions lead to vastly differ- 
ent dynamics of the CTQW. 

One way of quantifying the global efficiency of quantum 
walks is by the average probability of a walker to return to 
or stay at the origin [14]. In the classical case this quantity 
depends only on the eigenvalues, or more generally on the 
density of states (DOS), of the underlying system and not on 
its eigenvectors. Quantum mechanically, there exists a lower 
bound to the quantum mechanical return probability. This 
lower bound also depends only on the DOS ITil [Tsll . and, in 
most cases considered here, is a good measure of the transport 



efficiency, since the global temporal behaviour of the lower 
bound is similar to the one of the full expression for the return 
probability which also requires the eigenstates. 

Depending (mainly) on the topology of the system, the 
DOS shows very distinct features. There is a large variety of 
complex classical systems, ranging from glasses to proteins, 
showing anomalous transport 1 16], which can be related to the 
structure of the DOS. For instance, CTRW over small- world 
networks have been shown to be super-diffusive ITtIi . As we 
proceed to show, for quantum walks especially the degenera- 
cies of the eigenvalues and the number of degenerate eigen- 
values determine the transport efficiency. 

II. QUANTUM WALKS ON NETWORKS 

We start by considering quantum mechanical transport pro- 
cesses on discrete networks, which are a collection of N con- 
nected nodes. A connectivity matrix A = (Aij) can be 
assigned to every network. The non-diagonal elements Aij 
equal — 1 if nodes i and j are connected by a bond and oth- 
erwise. The diagonal elements An equal the number of bonds, 
fi, which exit from node i. 

A. Transition probabilities 

Classically, a CTRW is governed by a master equation for 
the conditional probability, pkj (t), to find the walker at time 
t at node k when starting at node j.|18] The transfer matrix 
of the walk, T = (Tkj), is, in the simplest case where the 
transmission rates 7 of all bonds are take to be equal, related 
to the connectivity matrix by T = — 7A (we assume 7 = 1 in 
the following). 

CTQWs are obtained by identifying the Hamiltonian of the 
system with the classical transfer matrix, H = — T |[3l [T2i[T3ll . 
The states | j) associated with the nodes j of the network form 
a complete, ortho-normalised basis set of the whole accessible 
Hilbert space, i.e., {k\j) ~ Skj- A state \ j) evolves in time as 
\j{t)) = \J{t)\j), where U(t) — exp(— iHt) is the quantum 
mechanical time evolution operator (we have set m = 1 and 
h= 1). 

The transition amplitude ak.j (t) from state | j) at time to 
state |fc) at time t reads then ak,j{t) ~ {k\XJ{t)\j) and obeys 
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Schrodinger's equation. Denoting the eigenvalues of the 
Hamiltonian H — —T hy En (n — 1, . . . ,N) and the ortho- 



Quantum mechanically, the average is given by 



normalised eigenstates by \ipn), such that J2n \i^n){i^n 
the quantum mechanical transition probability is 



(A:|e-'^"*|V„)(^nb) 



1, 
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By using the Cauchy- S ch warz inequality we obtain a lower 
bound for 7f(i) ifTlflsll. 



B. Long time limit 
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Quantum mechanically the unitary time evolution prevents ^^^.j^ q therefore get 



TTfe from having a definite limit for t oo. In order to 
compare the classical long time probability with the quantum 
mechanical one, one usually uses the long time average (LTA) 

in 

1 

Xk,j = lim - / dtnk,j{t), (2) 
= Y.^E^^E^{k\i^n){Mj){M^){'^^\k), (3) 

where 5E„,Em = 1 if En = and 5e„,e,,^ = otherwise. 
Some eigenvalues of H might be degenerate, so that the sum 
in eq. (O can contain terms belonging to different eigenstates 
and iV'm). 

We can use the Cauchy^Schwarz inequality to obtain a 
lower bound for the LTA 11511 . such that the time integral in 
eq. © fulfills 



^ di |afe,, (t)|2 >!. dt 



This results in 

Xk,j > Jim 



T^oo 



(4) 



(5) 



The only term in the sum over n in eq. Q which survives 
after integration and taking the limit T ^ oo is the one with 
El = 0. The corresponding eigenvector can be written as 



l/ViVEf^ib) m. Since {k\i,i) - (ViIj) 



1/VN we get with eq. ^ 



xfcj>iWi)(V'iij>r = ^- 



(6) 



C. Averaged transition probabilities 

Quantum mechanically as well as classically, we can cal- 
culated properties from the probability distributions, which 
solely depend on the eigenvalues of H or T, respectively, but 
not on the eigenstates. Classically, there exists a simple ex- 
pression for the average return probability to the initially ex- 
cited node (see, e.g., ref. [i20il ): 
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7f{t)> 
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In analogy to the classical case, the lower bound de- 
pends only on the eigenvalues and not on the eigenstates of 
H. 



III. QUANTUM WALK EFFICIENCIES 

Equations ^ and ( fTOb can be used to quantify the efficiency 
of the two transport processes [ff?]: In the classical case a 
quick decrease of must result - on average - in a quick 
increase of the probability for the walker to be at any other 
but the initial node. Thus, the transport away from the initial 
node 2 is the more efficient the quicker the decrease of vit). 
For instance, it has been shown for small-world networks that 
follows a stretched exponential rather than a power law 
as for regular networks, which gives rise to a quicker decrease 
of and thus to a super-diffusive behaviour |17]. Quan- 
tum mechanically, 7f(t), as well as the lower bound 
will show strong oscillations due to the unitary time evolution. 
However, one can use the envelope of this oscillations as a 
measure for the efficiency, for which the same arguments as in 
the classical case apply. For regular rf-dimensional networks it 
is straightforward to show that the envelope of decays 
as i^'', whereas decays as t^''/^ 11 1411 . Since the decay 
of is much quicker than the one of the quantum 

walk is more efficient than the classical random walk. 

For finite systems, p(t), 7r(i), and |cS(t)P do not decay ad 
infinitum. Classically, the averaged probability will drop in 
the course of time to the equipartition value of since we 
always have the eigenvalues i?i = which is non-degenerate, 
i.e.. 



1 iv 



-E„t 



(11) 



Quantum mechanically, both, 7f(t) and |a(t)p, will be oscil- 
lating about a value given by the LTA of TT{t) or its lower 
bound respectively. The long time average of 7f(t) 

follows as 



lim 



dt Tf{t) 



-E 

TV ^ 

3,n 



(12) 
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which still depends on the eigenstates \ipn)- Using again the 
Cauchy-Schwarz inequality to obtain a lower bound for x we 
get with eq. ( fTol i 
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(13) 



These long time averages give indications on the overall per- 
formance of quantum walks: If the LTA are larger than the 
equipartition value of on average most of the probabil- 
ity will remain at the initial node. 

The argument of using the lower bound |a(t)P instead of 
7f{t) goes only one way: Only if the envelope of ex- 
ceeds the one of the classical probability p{t) at any(most) 
times is the quantum walk less efficient than the classical ran- 
dom walk. However, if the envelope of | a (i) p lies below 
in general, little can be said about 7f{t). It might also happen 
that the envelope W{t) lies above the classical curve, although 
in the examples considered below this is not the case. 




time t 



FIG. 1 : (Colour on-line) p{t) and | Q (t) P with the appropriate scaling 
f^^^^ and t~^, respectively, for a ring of size A'' = 1000. 



IV. DEGENERACffiS OF EIGENVALUES 

The eigenvalues of a large variety of networks will be 
degenerate. By denoting the degeneracy of En by D„ = 
D{En), we can recast eqs. (|7]i and (fTol i into 



nit) > 



-T Dn 
N ^ 



respectively. As we will show, there is a profound difference 
between the temporal behaviour of the two quantities depend- 
ing on the eigenvalues and their degeneracy. For classical 
transport processes, the long time behaviour is dominated by 
the smallest eigenvalues, no matter what their degeneracy is 
- large eigenvalues with high degeneracies will eventually be 
suppressed by the exponential. In the quantum case, however, 
the degeneracies become important due to the unitary time 
evolution. 



the sum are those for which £"„ = Em- Then, the sum will be 
of order ©(iV) and therefore \a{t)f will be of order 0(1 /TV). 

As an example, fig. [T] shows p{t) (solid black line) and 
|a(t)p (short dashed red line) for a regular network of iV = 
1000 nodes. For regular networks, almost all eigenvalues are 
two-fold degenerate, for odd(even) N there is one(two) non- 
degenerate eigenvalue. Note that the lower bound is exact in 
this case, i.e., = W{t) 111- 

While propagating without interference, the envelope of 
|a(t)|^ decays as t^^ (dotted blue line), which is faster than 
in the classical case, where p{t) ^ t^^^^ (long dashed green 
line) [14]. Therefore, the quantum walk can be considered 
to be more efficient than the corresponding classical random 
walk. After the interference sets in at about t ~ 500 — N/2, 

\a(t)\'^, (15) fluctuates about the LTA x = (2A^ - l)/N'^ Ix = 

(2A^ — 2)/N'^] for odd(even) N, which is indeed of the same 
order as the classical long time equipartition value 1 /N. How- 
ever, the classical plateau value is reached at a much later time 
at about t sa lO''. 



(14) 



B. One highly degenerate eigenvalue 

In contrast to the uniform degeneracy, we consider now 
one eigenvalue, Ei, whose degeneracy, Di, is of order 0{N) 
whereas the others are of order 0(1) or less. By writing 



A. Uniform degeneracy of eigenvalues 

We start with networks, whose eigenvalues have degenera- 
cies of the order 0(1). Then we rewrite the lower bound as 



a(t) 



1 

N 



J2 Dne-' 

E„^Ei 



E„t 



we obtain, up to order 0{1/N'^), 



(17) 



E„,E^>E„ 



D„Dm cos[{En - Em)t]. (16) 



Now, the product will be of order 0(1), therefore, 

all terms in the sums contribute to |S(t)p. The shape of the 
decay depends on the functional form of the degeneracy. If we 
assume _D„ ^ E^" , we obtain p{t) ^ t~'^, whereas |a(t)|'^ ^ 
t~^^ Iil4i1 . At large times t the only terms which contribute to 



E 

E,,^Ei 



Dn2cos[{Ei- En)t] }. (18) 



The first term on the right-hand side of eq. ( fTSl ) is of order 
0(1) and the second term [with the same argument as for 
eq. ( fTSb l is of order 0(l/iV). Therefore, for few highly de- 
generate eigenvalues, the lower bound \a{t)\ will not show 
a decay to values which fluctuate about 1/A^ but rather will 
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fluctuate for all times about 1 — Also 7r(t) will not de- 

cay but fluctuate roughly about the same value since \a{t) \ is 
a lower bound. 
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FIG. 2: (Colour on-line), pit), 7v{t), and for a star with 

A'' = 51 nodes. The inset shows a close-up of 7r(t) and 
in the same time interval. 



Take as a simple example a star-shaped network, having one 
core node and — 1 nodes directly connected to the core but 
not to each other The eigenvalue spectrum has a very simple 
structure, there are 3 distinct eigenvalues, namely Ei — 0, 
E2 = 1, and E3 = N, having the degeneracies Di = 1, 
D2 ~ N — 2, and — 1, respectively. Therefore we get 



l + (Ar-2)e-'* + e"'(^-2)t 



(19) 



(20) 



1/N. 

for a star-shaped 



Obviously, only the term \{N ~ 2) exp(-it)p/Ar2 = (iV - 
2)^/N^ in eq. ^ is of order 0(1). All the other terms are of 
order 0{1/N) or 0{1/N'^) and, therefore, cause only small 
oscillations (fluctuating terms) about or negligible shifts (con 
stant terms) from {N - 2)^/N^ sa 1 

Figure |2] shows p{t), Tf{t), and | 
network with = 51 nodes. Here, not only (short 
dashed red line) but also Tf{t) (long dashed green line) fluc- 
tuate about a value close to one, see also the inset of fig. E] 
Classically, p{t) (solid black line) decays to the equipartition 
value 1/A^, although the exact form of the decay is different 
than for regular networks. Thus, the quantum walk has - on 
average - a large probability to return (or stay) at the initial 
site, and in this sense is less efficient than its classical coun- 
terpart. 



C. Two highly degenerate eigenvalues 

The next step is having two highly degenerate eigenvalues. 
El and £'.,„, with degeneracies Di and D„i of order 0{N/2). 



This changes the properties of . From 

ait) 



1 

N 



Di e-'^'* + Dme-'^-"' 



E 



Dr, e" 



(21) 



we get, up to order Oil/N), 
\a{t)? ~ -^[D^i+Dl^+2DiD,,, cos[{Ei-E„,)t]y (22) 

Here, the terms Df , D'^-^, and DiDm are of order 0{N'^ /A); 
the cosine fluctuates in the interval [0, 1]. Therefore, in con- 
trast to eq. ( fTSl l whose fluctuations are of the order 0{\/N), 
eq. (|22] | will also fluctuate in the full interval [0, 1]. As before, 
there is no envelope of |^ which decays with time. 




time t 

FIG. 3: (Colour on-line) p{t), Tv{t), and |a(t)|^ for a star with 50 
arms of length 2, i.e., N — 101 nodes. 

An example of such a system is a star with 2 nodes on 
each arm. There is one zero eigenvalue Ei = 0, two non- 
degenerate eigenvalues £^2, 3 = (iV+5±ViV^ -6N + 25)/4, 
and two eigenvalues £^4 5 = (3 ± \/5)/2, whose degeneracy 
^ D5 ^ {N - 3)'/2, i.e., it is of order 0{M), where 
M — {N — l)/2 is the number of arms. Then, for large N, 
where we assume N — 3 ^ N,eq. ( |22] | becomes 

1 



1 



cos 



{V5t) 



(23) 



Figure [3] shows the temporal behaviour of p{t), Tf{t), 
|a(t)p, and of eq. ( l23T l for a star with 50 arms of length 2, i.e., 
N = 101 nodes. While the behaviour of (black solid line) 
is similar to the one for the simple star (fig.|2]i, there are obvi- 
ous differences for 7f(t) (long dashed green line) and |a(i)p 
(short dashed red line). Although, as for the star-shaped net- 
work, neither 7f(t) nor show a decay to some value 
comparable to the classical equipartition value 1/iV, the os- 
cillations become larger for both quantities. While Tf{t) fluc- 
tuates roughly in the interval [0.2, 1], fluctuates in the 
full interval [0,1]. Moreover, considering only the two highly 
degenerate eigenvalues results in an excelent agreement of 
eq. (l23T l (dotted blue line) with the full expression for 
Thus also here the quantum walks have a large probability to 
return or stay at the origin, which results in a low transport 
efficiency compare to the classical random walk. 
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D. Complex structures with highly degenerate eigenvalues 

So far we have considered simple systems. As we now turn 
to more complex structure, it will become evident that a lot 
can be related to those simple models. In particular, there are 
also large complex systems of eminent interest which have 
only few degenerate eigenvalues. 

One group of systems of interest are star-like hyper- 
branched structures. An example is the dendrimer, which 
has been experimentally realised as one macromolecule and 
which has interesting applications, such as drug delivery and 
as a light-harvesting antenna lEsIl . The nodes of the dendrimer 
can be grouped into generations g, which are concentric about 
the core node. In our case the core {g = 0) has 3 emanating 
bonds connecting three additional nodes to it (g = 1); every 
node in generation g > 1 is connected by one bond to two 
other node in generation g + 1. Therefore, total number of 
nodes grows exponentially with the total number of genera- 
tions G, i.e., N = 3 • 2*^ — 2, whereas the number in one 
generation g is given by A^g = 3 • 2^^^. 

The eigenvalues of dendrimers can be calculated recur- 
sively 1 24], which yields one eigenvalue Ei — and G + 1 
non-degenerate eigenvalues. Additionally, there are N—G—2 
degenerate eigenvalues, whose degeneracies increase with the 
number of generations. In the first generation g, where the 
eigenvalue appears, it is two-fold degenerate, while for larger 
generation the degeneracy grows as 3 • 2^^^^^. The first 
appearing degenerate eigenvalues are £'2 = 1 and £"34 = 

2 ± V3, whose degeneracies are D2 — D3 = D4 = 3 • 2'^~^ 
for G > 2 and 3 • 2^^^ for G > 3. Therefore, roughly 

3 • 3 • 2'^^^ « A^/3 eigenvalues, i.e., 1/3 of the total num- 
ber of eigenvalues, are given by E2 and £3,4. 




FIG. 4: (Colour on-line) p{t), \a{t)\'^, and the approximate value for 
|a(t)P [eq. l l24b 1 for a dendrimer of generation G = 10, i.e., with 
N = 3070 nodes. 

Figure |4] shows p{t) (soHd black Une) and (short- 
dashed red line) for a dendrimer of generation G = 10, i.e., 
with N ~ 3070 nodes. Classically, p{t) decays in a similar 
fashion as in figs. |2] and |3] In the quantum case, the lower 
bound |S(t)p strongly oscillates about a value which is close 
to ~ 0.2. We stress again that |a(i)p is the lower 

bound to 7f (t) and especially reproduces the maxima of Tf{t) 



quite well lUSll . Thus, also 7f{t) > will be restricted to 

even higher values, see also fig. 8 in ref. Ill5i1 . 

When restricting ourselves to only the three most highly 
degenerate eigenvalues, the dendrimer resembles a star-like 
structures. From the eigenvalues we obtain (for G > 3) 



\a{t)\^ « -^<^ l + 4cos(\/3t) 



cos 



(Vst) + cos(t) 



>, (24) 



where Af is an appropriate normalization constant. Figure |4] 
also shows the temporal behaviour of eq. (|24] |. see long- 
dashed blue line, which is very similar to the full expression 
for |a(i)p. Hence, the general behaviour is indeed dominated 
by the most degenerate eigenvalues, there are only slight de- 
viations due to the remaining ones. From |a(t)p as well as 
from eq. ( l24l l. we find also here that quantum walks over den- 
drimers are much less efficient than their classical counterpart. 

There is certainly a large variety of networks, on which 
the quantum mechanical transport behaviour remains to be 
investigated. Prime examples are fractals. Classically it has 
been shown that the return probability depends on the spectral 
(fracton) dimension of the fractal, see, e.g., ll25ll . The eigen- 
values of certain fractals can also be calculated recursively 
1^6] . In some cases this gives rise to highly degenerate eigen- 
values. 



V. CONCLUSIONS 

We have analyzed the efficiencies of quantum walks over 
different discrete networks by relating the density of states to 
the average return probability of the walk. Depending on the 
connectivity of the networks, the spectra of the Hamiltonians 
are vastly differerent. The difference between classical ran- 
dom and quantum walks over networks can be summarised 
as follows: In the classical case the eigenvalues of the trans- 
fer matrix itself are governing the average return probabilities 
p{t). Quantum mechanically, the degeneracies of the eigen- 
values of the Hamiltonian dominate the temporal behaviour 
and the long time averages of the average return probability 
7f{t) and its lower bound After giving general ar- 

guments corroborating this statements, we illustrated this by 
simple examples, like regular and star-like networks and ex- 
tended those to more complex structures, like Cayley trees or 
dendrimers, which also have few highly degenerate eigenval- 
ues. 
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